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Abstract 

We describe a search for the pair production of first-generation scalar and vector leptoquarks in 
the eejj and evjj channels by the D0 Collaboration. The data arc from the 1992-1996 pp run at 
y/s = 1.8 TcV at the Fermilab Tevatron collider. We find no evidence for leptoquark production; 
in addition, no kinematically interesting events are observed using relaxed selection criteria. The 
results from the eejj and evjj channels are combined with those from a previous D0 analysis of the 
vfjj channel to obtain 95% confidence level (C.L.) upper limits on the leptoquark pair-production 
cross section as a function of mass and of /3, the branching fraction to a charged lepton. These 
limits are compared to next-to-leading-order theory to set 95% C.L. lower limits on the mass of 
a first-generation scalar leptoquark of 225, 204, and 79 CeV/c^ for j3 = 1, ^, and 0, respectively. 
For vector leptoquarks with gauge (Yang-Mills) couplings, 95% C.L. lower limits of 345, 337, and 
206 GeV /(? are set on the mass for (3=1,^, and 0, respectively. Mass limits for vector leptoquarks 
are also set for anomalous vector couplings. 



4 



I. INTRODUCTION 



A. Leptoquarks 

Leptoquarks (LQ) are exotic particles that couple to both leptons and quarks and carry 
color, fractional electric charge, and both lepton and baryon numbers |^ . Although the pat- 
tern of three generations of doublets of quarks and leptons suggests leptoquarks as a possible 
reason for an underlying unity, they are not required in the standard model. Leptoquarks, 
however, do appear in composite models, technicolor theories, grand unified theories, and 
superstring-inspired Eg models. They are not part of the minimal supersymmetric (SUSY) 
standard model, but can be accommodated in certain extended SUSY models. Leptoquarks 
can be scalar (spin 0) or vector (spin 1) particles. In many models, both baryon and lepton 
numbers are conserved, allowing low-mass leptoquarks to exist without mediating proton 
decay. 

Leptoquarks with universal couplings to all flavors would give rise to flavor- changing neu- 
tral currents and are severely constrained by low-energy experiments. We therefore assume 
in our analysis that there is no intergenerational mixing and that, e.g., first-generation lep- 
toquarks couple only to e or v^. and to m or (i quarks. In most models containing leptoquarks, 
each leptoquark species has a fixed branching fraction to i'^q: /3 = 1, | or 0. Models with 
intergenerational mixing or extra fermions can have any value of /? between and 1. 

The HI and ZEUS experiments at the e^p collider HERA at DESY have published lower 
limits on the mass of a first-generation leptoquark that depend on the unknown leptoquark- 
lepton-quark coupling, A [Q, ^, ||, |^, ^ |^, ||, |, Pair production of leptoquarks, nearly 
independent of the value of A, could occur in e^e~ collisions via a virtual 7 or Z in the 
s-channel and in pp collisions via an intermediary gluon. Experiments at the CERN e^e~ 
LEP collider |TT, |T^, |TB|, H] and at the Fermilab Tevatron |T^, [T^, have searched for 



leptoquark pair production and have set lower limits on the masses of leptoquarks. 

In February 1997, the HI and ZEUS experiments reported an excess of events at high 
|jl8|, |T^. A possible interpretation of these events is the resonant production of first- 
generation leptoquarks at a mass (Mlq) near 200 GeV/c^ pO|. Additional data collected in 
1997 did not confirm this excess (For a recent review of leptoquark phenomenology 

and the status of leptoquark searches at HERA and the Tevatron, see Ref. |]l|.) 

B. Leptoquark Production at the Tevatron 

At the Tevatron, pair production of leptoquarks can proceed through quark- ant iquark 
annihilation (dominant for Mlq > 100 GeV/c^) and through gluon fusion, and is therefore 
independent of the LQ-e-g Yukawa coupling A. Pair production of first-generation lepto- 
quarks can result in three final states: two electrons and two jets (eejj); one electron, a 
neutrino, and two jets (eujj); or two neutrinos and two jets {vvjj)- The decay branching 
fractions in the eejj, ez/jj, and wji channels are 2/?(l — /?), and (1 — respec- 
tively. The cross section for pp LQ LQ eejj is therefore proportional to We use 
the next-to-leading-order (NLO) calculation of the pair-production cross section of scalar 
leptoquarks to compare our experimental results with theory. This calculation has a the- 



oretical uncertainty of about 15% which corresponds to the variation of the renormalization 
scale /i used in the calculations between /i = 2Mlq and /i = |Mlq. For vector lepto- 
quarks, NLO calculations are not yet available, and we therefore use the leading-order (LO) 
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pair-production cross section |2^. We consider three gluon couplings: Yang-Mills gauge 
couplings {kg = = 0), minimal vector anomalous couplings {kq = 1 and Ag = 0), and 
the anomalous couplings that yield the minimum cross section for 150 GeV/c^ leptoquarks 
at = 1.8 TeV {kg = 1.3 and Ag = -0.21) |3[. 



II. D0 DETECTOR AND TRIGGERING 

The D0 detector is a general-purpose detector consisting of three major systems: a central 
tracking system, a uranium/liquid- argon calorimeter, and a muon spectrometer. These are 



described in Ref. [Q. The features most relevant to this analysis are summarized below. 

The central tracking system has a cylindrical vertex drift chamber, a transition-radiation 
detector, a cylindrical central drift chamber, and drift chambers in the forward regions. 
The tracking system is used to determine the longitudinal {z) position of the pp interaction 
and to find tracks associated with electrons and muons. Information from the transition- 
radiation detector helps separate electrons from charged pions. The calorimeter consists of 
a central calorimeter (CC) that covers the detector pseudorapidity region Ir^dctl < 1-2 
and two end calorimeters (EC) that cover 1.5 < Ir^dctl < 4.2. Scintillation counters located 
in the intercryostat region provide information about jets for 1.2 < |?7det| < 1-5. The 
electromagnetic (EM) and hadronic calorimeters are segmented into cells in pseudorapidity 
and azimuthal angle (0) of size A?7det x = 0.1 x 0.1 (0.05 x 0.05 at EM shower maximum). 

The Main Ring synchrotron lies above the Tevatron beam line and passes through the 
outer section of the central calorimeter. Protons used for antiproton production pass through 
the Main Ring while the Tevatron is operating. Interactions in the Main Ring can cause 
spurious energy deposits in the calorimeter leading to false missing transverse energy {]^t) 
in collected events. Certain triggers are rejected when the protons are being injected into 
the Main Ring, every time the Main Ring beam passes through the detector, and during the 
subsequent "calorimeter recovery" period; other triggers are rejected during injection and 
when the proton bunch is present, but accepted during calorimeter recovery periods (called 
a "minimal" Main Ring veto). Since all events are tagged with the state of the Main Ring 
at the time of collection, this rejection can be performed offline for triggers relying on less 
restrictive Main Ring requirements. 

D0 employs a three-level trigger system. Level uses scintillation counters near the 
beam pipe to detect an inelastic collision; Level 1 sums the EM energy in calorimeter towers 
of size A?7det x = 0.2 x 0.2. Level 2 is a software trigger that forms clusters of calorimeter 
cells and applies preliminary requirements on the shower shape. Certain triggers also require 
energy clusters to be isolated. 



III. EVENT RECONSTRUCTION AND PARTICLE IDENTIFICATION 



The D0 reconstruction program, D0RECO, processes the triggered data into events with 
kinematic quantities and particle identification. This includes finding interaction vertices, 
tracks, and jets, and identifying electrons and muons, each with loose quality criteria to 
reject poorly-measured objects. Additional requirements are then applied for each analysis. 
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A. Electron Identification 



Electron identification for tlie eejj and eujj analyses is very similar. Electron candidates 
are first identified by finding isolated clusters of energy in the EM calorimeter. These EM 
clusters are required to be in the fiducial volume of the detector, i.e., Ir^detl < 1-1 (CC) 
or 1.5 < |?7det| < 2.5 (EC). EM clusters with a matching track from the primary vertex 
are called electrons; those without a matching track are called trackless electrons. A track 
and an EM cluster in the CC match if the distance between the track and the EM cluster 
centroid is small: 



Ctrk 



' + — <io. 



where is the azimuthal mismatch, Az is the mismatch along the beam direction, and 
5x is the resolution for the observable x. In the EC, Az is replaced by Ar, the mismatch 
transverse to the beam. 

For the eejj analysis, at least one of the two electrons in an event is required to have 
a matching track. An electron track can be improperly reconstructed due to inefficiencies 
in the central tracking chambers or because of poor track/EM cluster matching caused by 
incorrect vertex information. Using trackless electrons restores some of this lost efficiency, 
but at the expense of increased background. They are not used in the evjj analysis. 

For electron candidates with a matching track, we apply a likelihood test based on the 
following five variables: 

• Agreement between the observed shower shape and that expected for an electromag- 
netic shower. This is computed using a 41-variable covariance matrix for energy de- 
position in the cells of the electromagnetic calorimeter (if-matrix @)- 



• The ratio of the shower energy found in the EM calorimeter to the total shower energy, 
the electromagnetic energy fraction (EMF), is required to be that expected for an EM 
shower. 

• A small track match significance, cxtrk, is required. 

• The ionization dE/dx along the track is required to be that for a single minimum- 
ionizing particle. 

• A variable characterizing the energy deposited in the transition-radiation detector is 
required to be consistent with the expectation for an electron 

To a good approximation, these five quantities are independent of each other for electron 
showers. For EM objects without a matching track, an if-matrix < 100 is required. 

All EM objects are required to have deposited most of their energy in the EM calorimeter 
(EMF > 0.9). We also require EM objects to be isolated, using the variable: 

^ _ Etotjn = 0.4) - EemCR- = 0.2) 
EEM(7^ = 0.2) 

where EtotiTZ = 0.4) and Eem(JZ = 0.2) are the total and EM energies in a cone of radius 
TZ = J (A?7)2 + (A0)2 = 0.4 or 0.2 centered on the EM cluster, where the pseudorapidity 



is measured with respect to the interaction vertex [^. For electrons with matching tracks, 
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we require T < 0.15. To reduce the multijet background by about 50% in dielectron data 
in which one electron does not have a matching track, we require that electron to have 
X < 0.10. The electron identification criteria are summarized in Table |. 

The electron Et resolution is a{ET)/ET = 0.01570(0.072 GeV^^VVte) ©0.66 GeV/Er, 
where © denotes a sum in quadrature. The resolution in rj and for an electron is excellent, 
less than 10" 



B. Jet Reconstruction 



Jet reconstruction p8[ is based on energy deposition in calorimeter towers (the calorimeter 
cells within Arj x = 0.1 x 0.1) with Et > 1 GeV. Starting with the highest-ii^r tower, 
the energy deposited in a cone of radius 7^ = 0.7 around the center of the tower is summed 
and a new energy-weighted center is determined. This procedure is repeated, using the new 
center, until the jet's direction is stable. Only jets with Ej- > 8 GeV are retained. The final 
direction of a jet is given by: 



6'jct = tan ^ 



/{j:.E^)^ + {j:.Eiy 



TE^ 




where the polar angle Q is measured relative to the interaction vertex, Er^ = Ei sin(^^j) cos(0i), 
Ey = Ei sin(^j) sin(0i), E^ = Ei cos(^i), and i corresponds to all cells that are within 71 = 0.7. 
Jets are required to have |?7dct| < 2.5 and EMF < 0.95. 

The measured jet energy is corrected for effects due to the underlying event and out-of- 
cone showering in the calorimeter. The transverse energy resolution for central jets (|?7det| < 
0.5) varies from a{ET)/ET = 0.154 for Et ^ 36 GeV to a{ET)/ET = 0.050 for Et ^ 
The resolution in both t] and for 50 GeV jets varies from approximately 
< 0.5 to approximately 0.06 for 2.0 < jr^detl < 2.5 and improves as the jet 



300 GeV p8 
0.02 for |r/det| 
energy increases. 

We use jets reconstructed with the large 71 = 0.7 cone size to decrease the number of 
final-state-radiation jets that are reconstructed separately from the parent jet and to improve 
the jet-energy and mass resolutions. Jets are ordered in descending value of Et, with ji, 
the leading jet, having the highest Et- 



C. Missing Transverse Energy 

Missing transverse energy is calculated as 
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TABLE I: Electron identification requirements. 



Requirement Electrons with tracks Electrons without tracks 

Fiducial volume l^ydetl < 1-1 or 1.5 < |?7det| < 2.5 |?7det| < 1-1 or 1.5 < |r?det| < 2.5 

Track match significance (Ttrk < 10 

Electromagnetic fraction EMF > 0.9 EMF > 0.9 

EM cluster isolation T<0.15 T<0.10 

EM cluster shape ff-matrix ^ < 100 

5-variable likelihood < 1.0 



where 

= -J2E,smie,)cosi<l),)-Y.AEi 

i j 
i j 

The first sum is over all cells in the calorimeter and intercryostat detector above the noise 
threshold, and the second is over the corrections in Et applied to all electrons and jets in 
the event. The J^^ resolution is approximately 4 GeV per transverse component [29| and 
grows as the amount of calorimeter activity increases. 



D. Vertex Finding 

The standard D0 vertex-finding algorithm uses tracks found in the central tracking sys- 
tem to locate the intersection of groups of tracks along the beam line. The group with the 
largest number of tracks is chosen as the primary vertex. However, since there is an average 
of 1.5 interactions per beam crossing, the hard-scattering vertex is not always chosen cor- 
rectly by this algorithm. Using the electron to verify or recalculate the vertex significantly 
improves this efficiency The electron revertexing algorithm uses the track that best 

matches an EM calorimeter cluster and then recalculates the position of the vertex based 
on this track. The z position of the vertex is calculated by fitting a straight line through 
the centroids of the EM cluster and the matching track. We require every event to contain 
at least one EM object with a matching track usable for revertexing. If both EM clusters 
have a matching track, the primary vertex is calculated based on information from both of 
them. The kinematic properties of the objects (electrons, jets, J^t) in the event, such as 
transverse energy and pseudorapidity, are then recalculated based on the new vertex. All 
further analysis is done using the recalculated quantities. 

Figure |l| illustrates the improvement in the resolution of the Z-boson mass, as well as the 
reduction in background due to vertex misidentification for ee) + 2j events, after the 
revertexing. Events in this plot are allowed to have one EM cluster without an associated 
track. 



IV. SEARCH STRATEGIES AND OPTIMIZATION 

The choice of variables, and the selection of their optimal values, for improving the ratio 
of signal to background events is at the heart of searches for new particles. We use two 
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FIG. 1: ee) + 2j data before (solid) and after (dashed) revertexing: (a) has a Hnear scale 

and illustrates the improvement in the Z-boson mass resolution after the revertexing; (b) has a 
logarithmic scale and shows the suppression of the background from vertex misidentification in the 
tails of the Z-boson peak. 



optimization techniques to aid in this selection: the random grid search method, which has 



been used by D0 in the measurement of the top-quark pair-production cross section pi 
and in the search for the supersymmetric partner of the top quark ||3^, and neural network 
analysis, which has been used by D0 in the measurement of the top-quark mass |^7|, ^ and 
in the determination of the tt-to-all-jets cross section [Bl, ISH . 



A. Additional Variables 



In addition to kinematic variables such as the transverse energies of electrons and jets and 
the used in standard analyses, we study other variables to determine their efficiency in 
separating signal from background. These include the energy sums, event-shape variables, 
invariant-mass variables, and mass-difference variables listed below. 

• Energy and transverse energy sums 





— sum of the Et 


of the two leptons 




(two electrons. 


or electron and neutrino {$t)) 




— sum of the Et 


of all jets 




— sum of the Et 


of the two leading jets 


Hi''' 


— sum of the Et 


of the three leading jets 


Sx 


= H^ + Hi 




q12 


= H^ + H^^' 




s 


— total energy in 


the event 
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Event-shape variables 

centrality (St/S) 
aplanarity of jets and leptons ^T\, ^ 



sphericity |36 



rms of the -Er-weighted distribution in jet rj pl | 
Invariant-mass variables 

Mee — dielectron invariant mass 

Mej — invariant mass of various electron and jet 

combinations 
M^'^ — electron-neutrino transverse mass 

Mass-difference variables for the eejj analysis 



-(Mlq) 



(5M,_ , J(Mlq, - Mlq)2 + (Mlq, - Mlq) 



M ^ Mlq 

m ^ Mlq, - Mlq, 
M (Mlq,+MlqJ/2 

6M _ Mlq, -Mlq, 



M ^(Mlq, + Mlq,)/2 

where Mlq, and Mlq, are the electron-jet invariant-mass combinations that are closest 
to each other, and Mlq is the hypothesized leptoquark mass. 

Mass-difference variable for the eujj analysis 
SM 



M 



(Mlq) 



mm 



Meji - Mlq 




Me,2 - Mlq 




Mlq ' 


Mlq ; 



where Meji and Mej2 are the invariant masses of the electron with the first and the 
second jet, respectively, and Mlq is the hypothesized leptoquark mass. 

Over fifty combinations of these variables were used in the random grid search and neural 
network studies described below to determine the optimal set of variables and selection 
criteria for the eejj and eujj channels. 



B. Optimization Criterion 

If first-generation leptoquarks with a mass of approximately 200 GeV/c^ exist, we want 
to achieve the highest-possible discovery significance. If there is no evidence of leptoquark 
production, we want to set the lowest possible 95% C.L. limit on their production cross 
section. Based on the Monte Carlo (MC) simulations of the signal and the background 
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estimates described below, we pursue a fixed-background strategy for our searcli. We opti- 
mize our selection criteria by maximizing the signal efficiency for 0.4 expected background 
events. This method leads to excellent discovery potential and a 67% probability that no 
background events will be observed. If no events are observed, the experimental limit has 
the advantage of being independent of the predicted number of background events and its 
uncertainty. 



Random Grid Search 

The random grid search method, which was implemented as the computer program 
RGSEARCH [0, helps determine the set of cuts that optimally separates signal from back- 



ground. In a standard grid search, the signal and background acceptances for some cutoff 
(a^cut) on a variable x are determined for all values between some minimum and maximum, 
Xmin and Xmax, respectively. A refinement of this technique is to use the MC signal to define 
the range of Xcnt- For each MC event, Xcut is set to the generated value of x, and the accep- 
tances for signal and background are determined for that Xcut- While running RGSEARCH, 
the value of a cutoff on a variable can be fixed or allowed to vary in some range. Mini- 
mum and/or maximum values for Xcnt can be preset or, alternatively, any values that are 
allowed for signal can be used in the search. In general, the search is multidimensional, 
and many combinations of variables, both fixed and varying, are studied to find an optimal 
set of requirements to impose on the data. Trigger thresholds and other criteria used to 
define the initial data sample are also imposed in all RGSEARCH trials. One of the results 
of an RGSEARCH trial is a plot of the number of expected signal events versus the predicted 
number of background events, normalized to the luminosity of the data sample, including 
detection efficiencies. 



D. Neural Network Analysis 

We also use three-layer feed-forward neural networks (NN) |3^, |3^ in the search for 



leptoquarks. For each combination of n variables, a network is trained using MC signal events 
(5") and an appropriate mixture of background events (B) to yield an output discriminant 
Pnn near 1 for signal and for background. For a sufficiently large sample of training 
events, when the trained network is applied to the data, the discriminant output from the 
neural network is approximately -s(^^b(x)' where S{x) and B{x) are the n-variable signal 
and background densities. This defines contours of constant probability for signal versus 
background in the n-dimensional space that represent the optimal functions separating the 
signal from the background. The discriminant then becomes a single variable that can be 
used to optimize the analysis for any desired signal to background ratio. 



V. THE eejj CHANNEL 

The study of the eejj channel is particularly important because it is the only channel 
sensitive to leptoquarks with (3=1. It is also sensitive to leptoquarks with /3 < 1; however, 
since both leptoquarks have to decay in the charged-lepton mode, the cross section for 
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leptoquark pair production and subsequent decay into the eejj channel is suppressed by a 
factor of 

Independent of the scalar or vector nature of leptoquarks, the analyses are very similar. 
In particular the data sample and the final event selection are identical. We describe the 
scalar leptoquark analysis first, in detail, and then the vector leptoquark analysis. 



A. The Data 

1 Triggers 

Events with two electrons satisfying the online trigger requirements listed in Table |I^ are 
used as the starting sample for the dielectron data sample. The total integrated luminosity 
for these triggers is 123.0 ± 7.0 pb~^, which corresponds to sample of 9519 events. The 
average trigger efficiency for the data in this analysis is (99.5 ± 0.5)%. 



2 Event Selection for the Base Data Sample 

We require two electrons with > 20 GeV and at least two jets with Ej. > 15 GeV. As 
described in Sec. only one of the electrons is required to have a matching track. Events 
containing an electron close to a jet {ATZe < 0.7) are rejected. Events whose dielectron 
invariant mass lies inside the Z'-boson mass window, 82 < Mee < 100 GeV/c^, are also 
removed. After identification, fiducial, initial kinematic, and Mge requirements, 101 events 
remain. We call these events the base data sample. 



B. MC Signal Samples 



Leptoquark pair production in the eejj channel can be modeled as the production of a 
pair of identical strongly-interacting particles, each of which decays into an electron and a 
jet. Monte Carlo events simulating the pair production of scalar leptoquarks are generated 
using ISAJET pO] for leptoquark masses from 80 GeV/c^ to 250 GeV/c^. The ISAJET samples 
are used only for calculating acceptances; the NLO calculation of Ref. 
production cross section. 



22 is used for the 



C. Background Samples 

The primary backgrounds to the eejj final state are from e~^e~ ( "Drell-Yan" ) production 
with two or more jets, tt production, and multijet events in which two jets are misidentified 
as electrons. 



1 Drell-Yan Background 

Drell-Yan (DY) events are generated using ISAJET in four mass ranges: 20-60 GeV/c^, 
60-120 GeV/c^, 120-250 GeV/c^, and 250-500 GeV/c^. For calculating the background, 
the DY+2j cross section from ISAJET is normalized to the observed number of events in the 
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TABLE II: The Level 2 triggers used in the eejj analysis. The runs listed correspond to different 
periods during Run 1 of the Tevatron (1992-1996). The transverse energy of an EM cluster is 
denoted by E^^. The number of events is that in the initial data set. 

Run Trigger Requirements Integrated Luminosity Number of Events 

Run lA E:™''™^ > 10 GeV 14.7 pb"! 1131 

Run IB E^^^ > 20 GeV, isolated 97.8 pb-^ 7500 

Run IC Ef^i > 20 GeV, isolated 10.5 pb^^ 



Z-boson mass peak after imposition of the kinematic criteria described above. The scaling 
factor is 1.7 ± 0.1 and reflects the fact that ISAJET does not provide the NLO corrections 
("K-factor") to the LO DY production cross section. The uncertainty in this background is 
20%, dominated by the 15% uncertainty in the jet energy scale. We estimate that the base 
data sample contains 66.8 ± 13.4 DY events. 



2 tt Background 



The tt dileptons MC sample is produced using HERWIG for rrit = 170 GeV/c 



The events are representative of all ee, e/i and fifi final states, including those from r decay. 
The sample of 101,339 events corresponds to an integrated luminosity of about 270 fb~^. 
The D0 measurement [^| of the ti production cross section has an uncertainty of 35%. 
This, when combined with the 15% uncertainty in the jet energy scale, leads to an overall 
uncertainty of 38% in the predicted number of ti events. The base data sample is estimated 
to contain 1.8 ± 0.7 ti events. 



3 Photon Background 

Direct photon production is the main source of real photons (observed as EM objects 
without associated tracks) in the eejj final state; its contribution is small and is taken 
into account when the multijet background is estimated. Other sources of photons, such as 
Wj + 2j production, are negligible for high-ET photons. 



4 Multijet Background 

The multijet background is estimated using data collected with a trigger that required 
three jets with Ej, > 10 GeV at Level 2. This trigger was prescaled and had an integrated 
luminosity of 0.936 pb~^. Two sets of events are selected from this trigger. Events in the 3j 
sample are required to have at least two jets with Ej, > 15 GeV and at least one additional 
jet with E^ > 20 GeV. Events in the 2j+EM sample have an EM object with ^ > 20 GeV 
rather than a third jet. 

The probabilities for a jet to be misidentified as either an electron or trackless electron 
are determined by comparing the number of candidates with E^ > Eq that pass standard 
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quality cuts in the 2j+EM sample and the total number of jets with Ej< > Eq in the 3j 
sample. The threshold is varied from 20 to 50 GeV, and the probabilities are stable for 
a cut value above 25 GeV, i.e., above the jet trigger turn-on. The probabilities for a jet to 
be misidentified as an electron with a track and without a track are measured to be: 

/tra.k = (3.50 ± 0.35) X 10-^ /no track - (1-25 ± 0.13) X IQ'^ 

and, within the uncertainties, are independent of the Et and pseudorapidity of the electron. 
These values are cross-checked using the ratio of 3j+EM and 4j events. This method of 
determining the misidentification probability automatically accounts for the direct photon 
background that is a part of the general "multijet" background. 

We then apply these misidentification probabilities to the weighted number of 4j events 
in the 3j sample. The weight assigned to each event is the number of jet permutations that 
can be used to misidentify a pair of EM objects. The backgrounds in the two samples, two 
electrons or an electron and a trackless electron, are estimated by multiplying the weighted 
number of events by /track 2/track/no track, respectively. We assign an uncertainty of 15% to 
these values, which reflects the variation of the misidentification probabilities as a function 
of any difference between the CC and EC, as well as certain jet trigger turn-on effects. 
The number of misidentified multijet events in the base data sample is estimated to be 
24.3 ± 3.6 events. 

5 Total Background 

The total background estimate for the base data sample is 92. 8± 13.8 events, in agreement 
with the 101 events observed in the data. 

D. Electron Identification Efficiencies 

There are approximately 300 Z-boson events remaining in the initial data sample after 
all requirements except those on the dielectron mass and for electron identification. This is 
sufficient to estimate the identification efficiencies for CC-CC, CC-EC, and EC-EC electron 
combinations. 

We plot the dielectron mass spectrum without any electron identification requirements 
beyond EM object reconstruction and subtract the multijet and DY backgrounds using the 
standard "side-band" technique. We then apply the electron identification requirements, 
again subtracting the backgrounds using the same side-band technique. The ratio of the 
background-subtracted number of Z bosons with the identification requirements to that 
without the identification requirements gives the efficiency per event. The efficiency is 
(74 ±3)%, (66 ±4)%, and (68 ±9)% for CC-CC, CC-EC, and EC-EC electron combinations, 
respectively. 

To calculate the average efficiency for leptoquark events, we find the relative fractions 
of the CC-CC, CC-EC, and EC-EC topologies. These are the same, within the errors, 
for leptoquark masses of 180, 200, and 220 GeV/c^, and equal (83 ± 2)%, (16 ± 1)%, and 
(1.1±0.2)% for CC-CC, CC-EC, EC-EC combinations, respectively. These fractions and the 
electron identification efficiencies give an overall electron identification efficiency of (73±4)% 
for leptoquark masses between 180 and 220 GeV/c^. 
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E. Event Selection Optimization 



1 Random Grid Search 



Extensive testing of combinations of the variables described in Sec. p.V A| shows that the 



use of a single variable, the scalar sum of the transverse energies of all the objects in the 
event, St, is the most powerful. Figure shows the St distribution for the base data sample, 
the predicted background, and a sample of 200 GeV/c^ leptoquark MC events. All of the 
leptoquark MC samples and the DY, tt, and 2j+EM background samples are used in the 
random grid search. The leptoquark events are used to set the trial threshold values for the 
different parameters. The number of predicted background events is determined using the 
three background samples. Shown in Fig. ^ is the predicted number of signal events versus 
the expected number of background events for three different RGSEARCH trials, where the 
samples have been normalized to an integrated luminosity of 123 pb~^, and the detection 
efficiencies, as well as the kinematic acceptance for the RGSEARCH thresholds, have been 
included. In these trials, the Et thresholds of the two electrons and the two jets are fixed 
to those in the base data sample. The thresholds varied are those for St alone, for ^(200) 
alone, and for these two variables together. When combined, St and the mass-difference 
variable yield a higher signal efficiency for very low values of expected background (less 
than 0.3 events), but the result is comparable to the use of St alone when the expected 
background is approximately 0.4 events. For the same expected background, using just the 
mass-difference variable leads to a 10% reduction in the predicted number of signal events 
compared to that using just St- Requiring St > 350 GeV leads to approximately 0.4 
expected background events (see Sec. |1V Bj). The highest value of St seen in the data is 



312 GeV, therefore no events pass this requirement. 



2 Neural Network Analysis 

The analysis based on the random grid search uses the linear sum St = + Hj,. 
However, it is possible that a function other than a simple linear sum is the optimal way 
to combine the two variables. The simplest way to compute this function is with a two- 
dimensional neural network. For this approach, we use a neural network with two input 
nodes (corresponding to the variables and Hj,), three hidden nodes, and one output 
node. The network is trained using the 200 GeV/c^ leptoquark MC sample as signal (with 
a desired network output Pnn = 1) and the observed admixture of DY, tt, and multijet 
events as background (with desired "Dnn = 0). Figure ^ shows the distribution of Pnn for 
the background, the 200 GeV/c^ leptoquark MC events, and the data. The discrimination 
between signal and background is good. 

Each value of Pnn defines a contour of constant probability between signal and back- 
ground in the {H^, H^p) plane. The expected distributions in x = ( ) space for a 
200 GeV/c^ leptoquark signal, the background, and the data are shown in Fig. |^. The 
contours corresponding to Pnn = 0.5, 0.8, and 0.95 are also shown. 

Selecting events with Pnn > 0.95 yields approximately 0.4 background events. The 
highest value of 2^nn in the data is 0.92 and no events survive the selection. The efficiency 
for identifying 200 GeV/c^ leptoquark events using the NN analysis is nearly identical to 
the efficiency found using the St analysis. Since the two methods give essentially equivalent 
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FIG. 2: St distributions for background (solid line histogram), data (solid circles), and Mlq 
200 GeV/c^ MC events (open triangles) for the eejj analysis. 



> 
o 

3 



0.1 0.7 0.3 0.-1 0.5 

Number of Background Events 



FIG. 3: Predicted number of Mlq = 200 GeV/c^ events versus the predicted number of background 
events for three RGSEARCH runs. The upper dotted line shows the variation with St- The lower 
dotted line shows the variation with ^(200). The structure (gaps) arises from an increase in 
acceptance for DY events. The more dispersed set of dots shows the result when both St and 
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(200) are varied. The density of points is irrelevant. 
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FIG. 4: Comparison of Pnn distributions for the predicted background (solid line histogram), 
200 GeV/c^ leptoquark events (dashed line histogram), and the data (hatched histogram). 
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FIG. 5: versus Hj. for (a) the predicted background, (b) 200 GeV/c^ leptoquark events, and 
(c) the base data sample. The curved lines correspond to Pnn = 0.5, 0.8, and 0.95 (from left to 
right). The area of a displayed square is proportional to the number of events in that bin, with 
the total number of events normalized to 123 pb^^. 

results for the final experimental limits, we use the simpler St analysis based on the random 
grid search described in Sec. |V E| . 

F. Checks 

1 St Distribution 

The modeling of the St distribution for high-mass DY events is checked by studying 
and Hj, separately, using data and MC events in the Z'-boson mass region. The aver- 
age value of for high-mass DY events (which provide most of the DY background) is 
approximately 250 GeV, corresponding to an Hj^ of approximately 100 GeV for St = 350 
GeV. The distribution of Hj^ for high-mass DY events is expected to be similar to that for 
Z + 2j events. Figure ^ shows the Hj, distribution for Z + 2j MC and data. In the region 
corresponding to the St cutoff for high-mass DY events {Hj, f» 100 GeV), the agreement 
is good. Disagreement between the Z + 2j MC events and the data at higher values of St 
stems from the LO calculations used in the simulation and does not affect the results of this 
analysis. 

In addition, we fit the Hij, distribution of the data to a sum of the DY and multijet 
backgrounds (the expected ti background is smaller than the uncertainties in the fit and is 
neglected). Figure |^ shows the Hj, distribution for the data and the result of the fit for the 
two backgrounds. The fit yields 77.5 ± 15.9 DY events and 24.6 ± 13.9 misidentified multijet 
events, for a total of 102 ± 21 events, in agreement with the 101 events in the base data 
sample and with the direct determination of the two dominant background contributions. 
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FIG. 6: The Hj^ distribution for Z + 2j data (solid circles) and MC (open triangles) in the Zboson 
mass region. For high-mass DY events, St ~ 350 GeV corresponds to Hj, ~ 100 GeV. 
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FIG. 7: Fit of the Hj. distribution in the eejj data to the sum of the DY and multijet backgrounds. 



2 Mass Fitting 

To improve resolution, rather than simply calculating the invariant masses of the electron- 
jet pairs, we use a kinematic fitter to reconstruct the mass of two identical particles that 
decay to electron+jet. The D0 fitting package, KFIT, is based on the bubble-chamber fitting 
program SQUAW 



The fitter balances the two electrons and the two leading jets against any extra jets and 
unclustered energy in the event by minimizing a to find the best fit solution. The takes 
into account the object resolutions (see Sec. |ITTD as well as the kinematic constraints. Three 



constraints are used in the fit: momentum conservation in the x and y directions for electrons, 
jets and unclustered energy, and the equivalence in the mass of the two leptoquarks. 

In each event there are two ways to associate the electrons and two leading jets (eiji, 62^2 
and eij2,e2ji). Fits for both configurations are performed and the configuration with the 
lowest is retained. The mass distribution for the background is found using the MC 
samples for DY and ti events; the multijet sample is not large enough to parametrize a 
smooth line shape, so a jet is used to simulate an electron in the fit. 
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Figure || shows St as a function of the fitted mass for the background, the 200 GeV/c^ 
leptoquark MC sample, and the data, before the St > 350 GeV requirement. The back- 
ground is centered at low St and low fitted mass and does not resemble the leptoquark 
signal. The data most closely resemble the expected background. Figure |^ displays the 
one-dimensional distributions in fitted mass for the three samples before the St cut and 
with a reduced St > 250 GeV requirement. The data and the predicted background are in 
good agreement. 



3 Varying the St Threshold 

Table shows a comparison between the predicted number of events from each of the 
three background sources, the total background, and the number of events observed in the 
data as a function of St threshold. The agreement between the predicted background and 
the data is excellent. 



G. Signal Studies 

1 Systematic Uncertainties 

The systematic uncertainties in the signal acceptance are obtained by comparing the 
results for scalar leptoquark samples generated using ISAJET and PYTHIA with different 
structure functions and renormalization scales. The uncertainty from the jet energy scale 
is determined by varying the calorimeter response to jets by one standard deviation. The 
systematic error in the signal varies from 17% to 13% for leptoquark masses between 120 



and 250 GeV/c . The systematic uncertainties are summarized in Table IV. 



2 Signal Efficiency 

The signal-detection efficiencies are determined using simulated scalar leptoquark events 
that pass the selection requirements and are shown in Table |V|. The uncertainties in the 
efficiencies include uncertainties in trigger and particle identification, the jet energy scale, 
effects of gluon radiation and parton fragmentation in the modeling, and finite Monte Carlo 
statistics. The overall efficiency ranges from 1% to 38.5% for leptoquark masses between 
80 GeV/c^ and 250 GeV/c^. 



H. Results from the eejj Channel for Scalar Leptoquarks 

Based on our observation of no events after requiring St > 350 GeV, we set a 95% 
C.L. upper limit on the leptoquark pair-production cross section using a Bayesian approach 
43| with a flat prior distribution for the signal cross section. Limits for different leptoquark 



masses are summarized in Table 0. As indicated before, to compare our experimental results 
with theory, we use the NLO calculation of the production cross section |2^. This cross 
section is tabulated for a wide range of leptoquark masses and has the value of 0.1841q;o26 pb 
for a 200 GeV/c^ leptoquark. The theoretical uncertainty corresponds to variation of the 
renormalization scale n used in the calculation from 2Mi^q to |Mlq. To set a limit on 
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FIG. 8: St versus the fitted mass for (a) background, (b) 200 GeV/c^ leptoquarks, and (c) tlie 
base data sample. Tlie area of a displayed square is proportional to the number of events in the 
bin. 



the leptoquark mass, we compare the theoretical cross section for /i = 2AfLQ with our 
experimental limit, resulting in Mlq > 225 GeV/c^ for a scalar leptoquark with /3 = 1 and 
Mlq > 176 GeV/c^ for a scalar leptoquark with /? = |. Figure |10| shows the experimental 
limit as a function of scalar leptoquark mass along with the predicted cross sections for 
/5 = 1 and (5 = \. The CDF collaboration has set a lower limit of Mlq > 213 GeV/c^ 

m 

1 



/3 = 1. When our result is combined with the CDF limit, a Tevatron mass limit of 
Mlq > 242 GeV/c2 is obtained for /3 
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I. Vector Leptoquarks 

Vector leptoquark events were generated for leptoquark masses from 100 GeV/c^ to 425 
GeV/c^ using a version of PYTHIA modified to include vector leptoquarks with various 



couplings. The distributions of the kinematic variables for scalar and vector leptoquarks are 
sufficiently similar that the same event selection can be used for both analyses. 

The identification efficiencies for vector leptoquarks for the three couplings considered are 
identical within their uncertainties, as shown in Fig. |ll]. To reduce the statistical uncertainty 
from the MC, we use the average identification efficiency of the three sets of MC events to 
set a single experimental limit on the cross section. This limit is then compared with the 
appropriate prediction for each coupling. 

The cross sections for vector leptoquark production have been calculated only to LO for 
three gluon couplings [^. For the scalar leptoquark case, cross sections calculated at NLO 
with \i = 2Mlq are approximately equal to those calculated at LO with Q'^ = Mlq. We 
therefore compare our cross section limit with LO calculations of vector leptoquark cross 
sections for this choice of scale. 
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FIG. 9: Distributions of fitted mass for the events in the base data sample (sohd circles), expected 
background (solid line histogram), and 200 GeV/c^ leptoquarks (hatched histogram) with (a) no 
cut on St and (b) a reduced threshold of St > 250 GeV. 

TABLE III: Comparison of the number of events expected from background with the number 
observed for the eejj analysis as a function of the threshold on St- 
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15 


225 


4.88 


1.73 





84 


7.45 


± 


1.06 


9 


250 


2.64 


0.99 





59 


4.22 


± 


0.59 


8 


275 


1.35 


0.60 





39 


2.34 


± 


0.32 


5 


300 


0.75 


0.35 





25 


1.35 




0.19 


3 


325 


0.31 


0.23 





16 


0.70 




0.09 





350 


0.18 


0.16 





11 


0.44 




0.06 





375 


0.12 


0.11 





07 


0.30 




0.04 





400 


0.07 


0.08 





04 


0.20 




0.03 






TABLE IV: Systematic uncertainties in the signal for the eejj analysis. 

Source Uncertainty (%) 

Particle identification 5 

Smearing in the detector 3 

Jet energy scale 11-2 (Mlq = 120-250 GeV/c^) 

Gluon radiation 7 

PDF and scale 7 

Monte Carlo statistics 2 

Luminosity 5 

Total 17-13 (Mlq = 120-250 GeV/c^) 
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TABLE V: Efficiency, background, 95% C.L. upper limit on the leptoquark pair production cross 
section (criimit)i and the NLO cross section (ctnlo) with fi = 2Mlq |22] for /3 = 1 as a function of 
leptoquark mass for the eejj channel. 



Mass 


Efficiency 


Background 


O'limit 


0"NLO 


(GeV/c2) 


(%) 


(Events) 


(Pb) 


(pb) 


80 


1.0 ±0.2 


0.44 ± 0.06 


2.9 


36.0 


100 


3.4 ±0.6 


0.44 ± 0.06 


0.80 


10.7 


120 


8.8 ± 1.4 


0.44 ± 0.06 


0.30 


3.81 


140 


14.4 ±2.1 


0.44 ± 0.06 


0.18 


1.54 


160 


20.9 ± 3.0 


0.44 ± 0.06 


0.13 


0.68 


180 


27.6 ± 3.8 


0.44 ± 0.06 


0.094 


0.32 


200 


33.2 ±4.0 


0.44 ± 0.06 


0.076 


0.16 


220 


36.1 ±4.4 


0.44 ± 0.06 


0.070 


0.080 


225 


37.7 ±4.5 


0.44 ± 0.06 


0.067 


0.068 


250 


38.5 ±4.7 


0.44 ± 0.06 


0.066 


0.030 



Figure 0(a) shows the experimental limits along with the three theoretical vector lepto- 
quark cross sections for the eejj channel for (3 = 1. Here, the experimental result yields a 
lower limit of Mlq > 340 GeV/c^ for the vector leptoquarks assuming Yang-Mills coupling, 
Mlq > 290 GeV/c^ for minimal vector coupling, and Mlq > 245 GeV/c^ for the coupling 
corresponding to the minimum cross section. Similarly, for /3 = | (Fig. |l2|(b)), our result 
provides a lower limit of 300 GeV/c^ for Yang-Mills coupling, 250 GeV/c^ for minimal vector 
coupling, and 210 GeV/c^ for the coupling corresponding to the minimum cross section. 

VI. THE eujj CHANNEL 

For < P < 1, leptoquark pairs can decay to eujj as well as to eejj. The eujj channel 
therefore allows us to extend the leptoquark mass limit to higher masses for < /3 < 1. Our 
optimization techniques for this analysis are similar to those we used for the eejj channel. 

As in the eejj channel, we use the same data sample for both the scalar and vector- 
leptoquark analyses. However, because the scalar-leptoquark analysis depends on a mass- 
based variable, and the vector leptoquark analysis is sensitive to higher masses than the 
scalar leptoquark analysis, the final event selection is slightly different. The scalar leptoquark 
analysis is described first, followed by the vector leptoquark analysis. 

A. The Data 

1 Triggers 

The data sample for this analysis corresponds to an integrated luminosity of 115 ±6 pb~^. 
Using events collected with the triggers shown in Table [VT| , the initial data sample contains 
95,383 events. 
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FIG. 10: Upper limit on the leptoquark pair-production cross section (triangles) from the eejj 
channel. The NLO calculations of Ref. ||2^ for /3 = 1 (upper band) and (3 = ^ (lower band) are 
also shown. The central lines correspond to /i = Mlq, and the lower and upper edges of the bands 
correspond to fi = 2Mlq and fj, = ^Mlq, respectively. 



2 Event Selection for the Base Data Sample 

We require one electron with a matching track with > 30 GeV, I^rp > 20 GeV, and at 
least two jets with > 20 GeV. Electrons with E^ > 20 GeV close to a jet {ATZe < 0.6), 
are "subtracted" from the jet in order not to double count the energy in the event. Since 
the J/^x threshold for this analysis is relatively high, we use a "minimal" Main Ring veto to 
increase the efficiency (see Sec. |VI D 2|) . 

To suppress the background from top-quark pair production, we apply a muon veto by 
requiring events to contain no well- reconstructed muons with p-r > 4 GeV/c |^6[. To reduce 
the multijet background when J^q- < 120 GeV, we require the vector to be isolated in </> 
from any jets (A0(j, J^q.) > 0.25). The effect of this requirement on a 180 GeV/c^ leptoquark 
MC sample and on the multijet background is shown in Fig. [1^. 



TABLE VI: The Level 2 triggers used in the ei^jj analysis. The transverse energy of an EM cluster 
is denoted by E^^. The number of events is that in the initial data set. 

Run Trigger Requirements Integrated Luminosity Number of Events 

Run lA E^^ > 20 GeV 11.2 pb^^ 9,862 

Run IB E^^ > 20 GeV, isolated 92.9 pb^^ 77,912 

> 15 GeV 

Run IC E^^ > 20 GeV, isolated 0.8 pb^^ 369 

> 15 GeV 

Run IC > 17 GeV, isolated 10.5 pb^^ 7,240 
^^^•^■^ > 10 GeV, > 14 GeV 
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FIG. 11: The efficiency for identifying vector leptoquarks for the three eouphngs in the eejj 
channel. The differences between the efficiencies are small relative to the uncertainties. 



After the above cuts, 1094 events remain in the data sample, primarily from W + 2j 
production. To remove these events, we require M^y > 110 GeV/c^, reducing our base data 
sample to 14 events. 



B. MC Signal Samples 

We use the ISAJET event generator followed by the full detector simulation via GEANT to 
model the leptoquark signal. Two thousand to five thousand events were generated in steps 
of 20 GeV/c^ for Mlq between 80 and 220 GcV/c^. We also use a PYTHIA MC sample at 
200 GeV/c^ for studying MG systematics and for cross checks. 



C. Background Samples 

As implied above, the dominant background to the eujj final state is W + 2j production. 
The other significant backgrounds are from tt production and multijet events in which a jet 
is misidentified as an electron and the energy is mismeasured, thereby introducing false 



1 tt Background 

The ti MG event sample contains all leptonic final states for rut = 170 GeV/c^. It was 
generated using HERWIG followed by GEANT detector simulation. The sample corresponds 
to an integrated luminosity of about 32 fb~^. 

Since top-quark events frequently contain muons from W ^ ixv and 6-quark decays, the 
muon-veto requirement provides an effective way to remove tt events. To determine the 
background due to top quark events, we apply all of the basic cuts except the muon and 
minimal Main Ring vetoes to the MG sample. 
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FIG. 12: The 95% C.L. upper limits on the vector leptoquark pair production cross section from 
the eejj channel and the LO predictions for Yang- Mills (YM) , minimal vector (MV) , and minimum 
cross section (MCS) couplings as a function of leptoquark mass for (a) /? = 1 and (b) /? = |. 



Because the reconstruction efficiency for muons in MC events is higher than that for real 
muons, GEANT overestimates the rejection factor against muons. The correction (between 
50% and 90%) to the efficiency depends on the run number (due to chamber aging and 
repair) and the pseudorapidity of the muon. After applying this factor and the efficiencies 
described below, we estimate that the data sample of 1094 events (before imposing the M^'^ 
cut) contains 12 ± 4 events. After requiring M^'^ > 110 GeV/c^, 2.0 ± 0.7 ti events are 
expected to remain in the base data sample of 14 events. 



2 Multijet Background 

The multijet background is estimated using the data samples and the misidentification 
probability of (3.50 ± 0.35) x lO"'' described in Sec. |V C4| . We select events from the 
multijet data sample that have at least three jets and > 30 GeV. To minimize luminosity 
dependence and the misidentification of primary interaction vertices, we use only those events 
that have a single interaction vertex within the fiducial region of the detector (|^;vtx| < 50 
cm) . To account for multiple interactions and multiple vertices, we apply a correction factor. 
The correction factor is determined by measuring the fraction of single-interaction events in 
the Z + 2j data sample as a function of luminosity, and then weighting this fraction with a 
luminosity profile of the multijet data stream. The correction factor is 2.2 ± 0.2. 

We next examine all three-jet combinations for each event. We treat each jet as an 
electron in turn and require each permutation to pass our electron and jet kinematic and 
fiducial requirements. Since the misidentification rate already accounts for the probability 
for a jet to be misidentified as an electron, we do not apply the electron identification 
criteria here. The multijet background is then defined by the product of the number of 
combinations that pass all criteria, the misidentification probability, and a factor that scales 
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FIG. 13: Effect of tlie requirement of acolinearity in $rp on (a) a 180 GeV/c^ MC leptoquark signal 
and (b) the multijet background. In (b), the dots show the distribution before imposition of the 
requirement; the open squares show the distribution after applying the MqV requirement. The 
acohnearity requirement is indicated by the sohd hues. 



the multijet sample luminosity to the luminosity of the data. There are 75 ± 15 events 
expected in the sample of 1094 events before the cut and 4.1 ± 0.9 multijet events after 
the > 110 GeV/c^ requirement. The uncertainty in the background accounts for the 
statistics of the multijet sample and for a 20% systematic error reflecting the variation of 
the misidentification probability with Et and pseudorapidity, as well as jet trigger turn-on 
effects and the uncertainty in the scaling factor. 



3 W + 2j Background 



For the W + 2j background, we use a sample of events generated with VECBOS ||46 
followed by ISAJET underlying-event modeling and GEANT detector simulation. This initial 
sample contains 227,726 events and corresponds to an integrated luminosity of approximately 
0.8 fb-i. 

For calculating the background, the number of MC W+2j events with M^" < 110 GeV/ 
is normalized to the observed number of events after subtracting the estimated ti and multijet 
backgrounds. A scaling factor of 0.22 ± 0.01 gives good agreement between the Monte 
Carlo and the data and is consistent with the value of 0.20 expected from cross section and 
efficiency calculations. 
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To check the normahzation, we repeat the comparison between the estimated background 
and the data for two additional thresholds on the ^5-: > 25 GeV and > 35 GeV. The 
agreement is again very good, showing that the fractional backgrounds are well-understood 
(the multijet background varies by a factor of six, from 115 to 20 events, between the two 
thresholds). 

The number of W + 2j events in the base data sample is estimated to be 11.7± 1.8 events. 
4 Total Background 

Figures |l5(a) and |14|(b) show the M^'^ and S]? distributions for the data sample and the 
background before the cut on M|>^. It is clear that we model the transverse mass distribution 
quite well up to 110 GeV/c^. The distribution is also well-described by the MC except 
for the small systematic offset of the prediction relative to the data. 

The total background estimate after basic requirements is 17.8 ±2.1 events, in agreement 
with the 14 events observed in the data. 

D. Efficiencies 

1 Trigger Efficiency 

Since events in the base data sample are required to have high electron Et and I^^, the 
trigger requirements listed in Table are very efficient. The EM part of the trigger has an 
efficiency of (99.5 ±0.5)%. 

2 Efficiency of the Minimal Main Ring Veto 

As discussed in Sec. |T|, additional Main Ring (MR) trigger requirements can be applied 
offiine to events collected using triggers with liberal MR requirements. For the eujj analysis, 
we apply a "minimal" MR veto to remove events that occurred during proton injection and 
when the proton bunch passed through the detector, while keeping events collected during 
the calorimeter recovery period. The efficiency of this veto is estimated using Z + 2j data 
collected using triggers with looser MR requirements than in the triggers used in the eujj 
analysis. First, the MR requirements for the eujj triggers are applied to the Z ± 2j data. 
The efficiency of the minimal MR veto is then calculated by comparing the number of events 
in the Z-boson mass peak before and after the additional minimal MR veto requirements 
are applied. The efficiency of this veto is (94 ± 1)% (i.e. 6% of the good events are removed 
along with a much larger percentage of background events). If the "calorimeter recovery" 
events were also removed, the efficiency would be reduced to about 90%. 

3 Muon-Veto Efficiency 

The efficiency of the muon veto is estimated using a sample of Z{^ ee)+2j events. Except 
for the additional electron, these events have a topology similar to that of leptoquark events 
in the eujj channel and should have a similar random muon track rate. The calculation is 
done using the number of events in the Z-boson mass peak before and after application of 
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FIG. 14: Comparison of the (a) M^'^ and (b) distributions for the evjj data (points with error 
bars) and the predicted background (soUd histogram) before imposing the requirement. 



the muon veto. Background under the Z boson is subtracted using the standard side-band 
technique. The muon-veto efficiency is (97 ± 1)%. 



4 Electron Identification Efficiencies 



Using the efficiencies described in Sec. [VD| for the eejj channel, the overall electron 
identification efficiency for leptoquark events in the ei^jj channel is (61 ±4)% in the CC 
and (54 ± 4)% in the EC. Since (93 ± 1)% of the electrons in the eujj final state are in the 
CC, the total electron identification efficiency, including tracking and quality requirements, 
is (60 ±3)%. 



E. Event Selection Optimization 



1 Random Grid Search 



We use a random grid search based on the Mlq = 180 GeV/c^ MC sample to select 
the optimal variables and thresholds for the eujj channel. Many different variables and 
combinations of variables (see Sec. [IV A\j were tested for their efficiency in retaining the 
signal and rejecting the background. The inputs to the RGSEARCH program are the MC 
signal samples and the W + 2j, ti, and multijet background samples described in Sec. |V1 C 



The combinations of variables that have the most discriminating power are then used in the 
neural network analysis. The most powerful variables for separating leptoquark signals from 
the background are and 



see Sec. WM - 
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2 Neural Network Analysis 



SM 
M 



We use a neural network with two input nodes (corresponding to the variables S]? and 
(Mlq)), five hidden nodes, and one output node. A separate network is trained for 
each MC signal sample (with a desired network output T'nn(^lq) = 1) and the expected 
admixture of W + 2j, tt, and multijet background events (with desired r'NN(^LQ) = 0). 
The expected rejection can be seen in Figs. ^ and |T6|. Figure ^ shows the two-dimensional 
distributions of ^(180) versus for the three individual backgrounds. Figures p!6|(a)-(c) 
show the same two-dimensional distributions for the total background, simulated leptoquark 
events with Mlq = 180 GeV/c^, and the data. The contours corresponding to constant 
values of X'nn(180) = 0.75, 0.85, and 0.95 demonstrate the level of separation achieved 
between the expected signal and the background. The distribution of X'nn(180) for the data 
is compared with the predicted distributions for background and signal in Fig. 0(d). The 
data can be described by background alone. The highest value of X'nn(180) observed in the 
base data sample is 0.79. 

Using the strategy described in Sec. [IV B| , we optimize the signal for a fixed background 
of approximately 0.4 events. In the low- mass range (Mlq < 120 GeV/c^), where leptoquark 
production rates are high, requiring S]? > 400 GeV is sufficient and leads to a background 
of 0.60 ±0.27 events, consistent with the desired background level. For Mlq > 120 GeV/c^, 
we use neural networks since they provide higher efficiency than an cut alone. For 
180 GeV/c^ leptoquarks, approximately 0.4 background events are expected for Pnn(180) > 
0.85. We choose the Pnn(Mlq) threshold to be a multiple of 0.05 rather than a value that 
yields exactly 0.4 background events; Pnn(180) > 0.85 corresponds to a background of 
0.29 ± 0.25 events. No events in the base data sample satisfy this criterion. Naturally, for 
leptoquark masses other than 180 GeV/c^, the requirement on T'nn(Mlq) is different. The 
expected background varies between 0.29 and 0.61 events and is listed in Table |VII| . No 
events from the base data sample pass any of these V^-^^Mi^q) thresholds. 

Rectangular cuts of S^^ > 350 GeV and ^(180) < 0.25 yield a total background of 0.4 
events. This also leaves no events in the data sample, but the signal efficiency is approxi- 
mately 10% lower for Mlq = 180 GeV/c^. 



Check 



As a check of our understanding of the background. Fig. 17 shows the distribution of 
Pnn(180) for the data and for the predicted background before the M^'^ cut. The agreement 
is acceptable. 



G. Signal Studies 

1 Systematic Uncertainties 

The systematic uncertainty in the signal efficiency varies from 25% to 8% for Mlq between 
80 and 220 GeV/c^. The sources and sizes of the systematic uncertainties are given in 
Table [Vlll . The uncertainties due to the jet energy scale and initial/final state radiation 



are significantly lower than in the eejj analysis due to the use of rather than the all- 
jets-based St discriminator. 
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FIG. 15: Distributions of ^(180) versus Sj? for the three individual backgrounds: (a) W + 2j 
events, (b) multijet events, and (c) ti events. The curves show neural net contours for Pnn(180) = 
0.75,0.85, and 0.95. 



2 Signal Efficiency 

The signal detection efficiencies are calculated using simulated leptoquark events that pass 
the selection requirements; they are shown in Table [VI 1| . The errors in the signal efficiencies 
include uncertainties in trigger and particle-identification efficiencies, the jet energy scale, 
effects of gluon radiation and parton fragmentation in the signal modeling, and finite MC 
statistics. 



H. Results from the evjj Channel for Scalar Leptoquarks 

We obtain a 95% C.L. upper limit on the scalar leptoquark pair-production cross section 
for /3 = I as a function of leptoquark mass. The results, based on a Bayesian analysis |^ 



are shown in Table |V14 The statistical and systematic uncertainties in the efficiency, the 
integrated luminosity, and the background estimation are included in the limit calculation, 
all with Gaussian priors. The 95% C.L. upper limits on the cross section for scalar leptoquark 
pair production in the eujj channel, corrected for the branching fraction of P = |, for 



various leptoquark masses are plotted in Fig. [T^ along with the NLO calculations p2| . 
The intersection of our limit with the lower edge of the theory band (renormalization scale 
/i = 2Mlq) is at 0.38 pb, leading to a 95% C.L. lower limit on the leptoquark mass of 
175 GeV/c2. 



I. Vector Leptoquarks 

As in the case of the eejj channel, vector leptoquark events were generated for Mlq 
between 100 GeV/c^ and 425 GeV/c^ using a version of PYTHIA modified to include vector 
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FIG. 16: Distributions of ^(180) versus S^j? for (a) the total background, (b) ten times the ex- 
pected signal from 180 GeV/c^ leptoquarks, and (c) the data, (d) The neural network discriminant 
for the signal (hatched histogram), the background (open histogram), and the data (points with 
error bars). The curves show neural net contours for Pnn = 0.95, 0.85, and 0.75. 



leptoquarks with different couplings. The distributions of the kinematic variables for scalar 
and vector leptoquarks are similar, and consequently, the same event selection is used for 
both analyses for Mlq < 220 GeV/c^. 

Neural networks for the evjj channel were trained on scalar leptoquark MC samples up 
to Mlq = 220 GeV/c^. Since vector leptoquark production cross sections are higher than 
scalar leptoquark cross sections, higher masses are of more interest. For vector leptoquarks 
with Mlq > 220 GeV/c^, we require S}j? > 400 GeV. This variable is one of the inputs to 
the neural network and provides good signal identification efficiency and a background of 
0.60 ± 0.27 events. 

Again, the identification efficiencies for vector leptoquarks for the three couplings agree 
within their uncertainties, as shown in Fig. [1^. Therefore, to reduce the statistical un- 




Dnn(180) 



FIG. 17: Comparison of the Pnn(180) distribution for the evjj data (points with error bars) and 
the predicted background (solid histogram) before the cut on M^. 
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TABLE VII: Efficiency, background, 95% C.L. upper limit on the leptoquark production cross 
section, and NLO cross section multiplied by the branching fraction with /j = 2Mlq [^2| for P = ^ 
as a function of leptoquark mass for the ei'jj channel. 



A/Toco 
iVicLSS 


Efficiency 


ijackground 


0" limit 


Zp[l — pjfTNLO 


(GeV/c^j 


(%) 


(Events) 


(Pb) 


(pb) 


80 


0.32 ±0.08 


0.60 ±0.27 


10.9 


18.0 


100 


1.15 ±0.21 


0.60 ± 0.27 


2.6 


5.34 


120 


2.45 ±0.33 


0.60 ± 0.27 


1.0 


1.90 


140 


6.65 ±0.96 


0.54 ±0.25 


0.43 


0.77 


160 


10.9 ±1.2 


0.61 ±0.27 


0.24 


0.34 


180 


14.7 ± 1.2 


0.29 ±0.25 


0.18 


0.16 


200 


19.4 ±1.7 


0.43 ± 0.27 


0.14 


0.08 


220 


21.5 ± 1.7 


0.41 ±0.27 


0.12 


0.04 



TABLE VIII: Systematic uncertainties in the signal for the eujj analysis. 

Source of Systematics Uncertainty (%) 

Particle identification 5 

Smearing in the detector 3 

Jet energy scale 10-2 (Mlq = 80-220 GeV/c^) 

Gluon radiation 4 

PDF and scale 5 

Monte Carlo statistics 25-3 (Mlq = 80-220 GeV/c^) 

Luminosity 5 

Total 25-8 (Mlq = 80-220 GeV/c^) 
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FIG. 18: The 95% C.L. upper limits on the cross sections for scalar leptoquark pair production 
from the eujj channel, and for all three channels combined, for /3 = ^, compared to the NLO 
prediction, clS Si function of leptoquark mass. 
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certainty in our analysis, we use the average identification efficiency of tlie tliree sets of 
MC events to set a single experimental limit on the cross section. As before, this limit is 
compared with the appropriate theoretical cross section for each coupling. 

Figure ^ shows the experimental limits along with the three theoretical LO vector lep- 
toquark cross sections [^] for the eujj channel for /3 = ^ and = MIq. For Yang-Mills 
coupling, the experimental lower limit on the vector leptoquark mass is 315 GeV/c^, for 
f5 = ^. For minimal vector coupling, the mass limit is 260 GeV/c^ for = i. For the 
coupling corresponding to the minimum cross section, the lower limit is 215 GeV/c^ for 

P = i 



VII. THE uujj CHANNEL 



To analyze the I'l'jj channel, we make use of our published search |^ for the supersym- 
metric partner of the top quark using just the 1992-1993 data sample. In that analysis, we 
searched for the pair production of top squarks that decay exclusively via a c quark and the 
lightest neutralino, ti — cXi, resulting in a final state with J^rp and two acolinear jets. 

Approximately 75% of the data were collected using a trigger whose primary requirement 
was I^rp > 35 GeV at Level 2; the balance had a threshold of 40 GeV. To ensure an un- 
ambiguous ]/]rp measurement, events were required to have only one primary vertex, reducing 
the sample to single interactions with an integrated luminosity equivalent to approximately 
7.4 ph-\ 

Events were required to have J^x > 40 GeV, two jets with Et > 30 GeV, and no isolated 
electrons or muons with Et > 10 GeV. In addition, the two leading jets were required to 
be acolinear (90° < A0(ji,j2) < 165°), and the 

was required not to be aligned with 
either the leading jet (10° < A0(ji,_^j.) < 125°) or the third or fourth leading jets (10° < 
A0(j3 4, ^j.)). Three events survived the selection criteria, consistent with the estimated 
background of 3.5 ± 1.2 events, primarily from W/Z+iets production. 

The efficiencies of the event selection for scalar leptoquarks with Mlq from 50 to 
200 GeV/c^ are calculated using signal MC events generated with the ISAJET generator 
and processed through the GEANT-based detector simulation. The systematic errors in the 
signal acceptance are calculated as in Ref. |^. The efficiencies, background, and cross 
section limits are shown in Table This analysis yields the limit Mlq > 79 GeV/c^ at 



the 95% C.L. for (3 = 0. 

The identification efficiency for vector leptoquark (generated using PYTHIA) and scalar 
leptoquark events with Mlq = 200 GeV/c^ are identical, within errors. Based on this 
comparison, and similar comparisons in the eejj and eujj channels, we use the experimental 
limit for scalar leptoquarks for vector leptoquarks in the i'l'jj channel. Comparison with the 
theoretical cross sections leads to 95% C.L. limits of Mlq > 206, 154, and 144 GeV/c^ for 
Yang-Mills, minimal vector, and minimum cross section couplings, respectively, for (3 = 0. 



VIII. GAP IN THE LIMIT FOR SCALAR LEPTOQUARKS 

In our analysis of the eejj and eujj channels, we use MC samples of leptoquarks 
with Mlq > 80 GeV/c^ but our analysis is optimized for leptoquarks with masses near 
200 GeV/c^. From the eujj analysis, we exclude /3 > 0.13 for Mlq = 80 GeV/c^. The mass 
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FIG. 19: The detection efficiency for vector leptoquarks in the ei'jj channel. 



limit from the z/z/jj channel for f3 = 0.13 is approximately 75 GeV/c^, leaving a small gap 
in our limit. 

To fill this gap, we examine further the 14 events in the base data sample in the eujj 
analysis. Making the very conservative assumption that all 14 events are due to leptoquark 
pair production, the 95% CL. upper limit on the cross section multiplied by the branching 
fraction and efficiency is 0.20 pb. This permits us to extend our exclusion region to include 
0.09 < P < 0.91 for Mlq = 80 GeV/c^ and 0.05 < /? < 0.95 for Mlq = 75 GeV/c^. To 
obtain the efficiency for Mlq = 75 GeV/c^, we scale the efficiency found for higher Mlq. 



TABLE IX: Efficiency, background, 95% C.L. upper limit on the leptoquark pair production cross 
section, and the NLO cross section with /i = 2Mlq [22] for /? = as a function of leptoquark mass 
for the i^i'jj channel. 



Mass 
(GeV/c2) 


Efficiency 

(%) 


Background 
(Events) 


Climit 

(pb) 


0"NLO 
(pb) 


50 




3.49 ±1.17 


328 


406 


60 


1.11 ±0.16 


3.49 ±1.17 


77.0 


162 


80 


9 1 c:+0.23 


3.49 ±1.17 


37.7 


36.0 


100 


3.90 ± 0.30 


3.49 ±1.17 
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3.49 ±1.17 
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FIG. 20: The 95% C.L. upper limits on the cross sections for vector leptoquark pair production 
from the evjj channel for /3 = i, and the LO predictions for the three couplings, as a function of 
leptoquark mass. 



IX. COMBINED RESULTS 



Combining |^ the limits from the eejj, ei^jj, and vyjj channels, we obtain 95% C.L. 
upper limits on the leptoquark pair-production cross section as a function of leptoquark 
mass and [3. The cross-section limits for [3 = \ are shown in Fig. |1^ for scalar leptoquarks 
and in Fig. ^ for vector leptoquarks. Table ^ lists the mass limits for (3 = 1, |, and 
for the types of leptoquarks studied. The lower limits on the mass of scalar leptoquarks as 
a function of /?, for all three channels combined, as well as for the individual channels, are 



shown in Fig. Figure ^ shows the exclusion contours from the individual channels and 
the combined result for vector leptoquarks with Yang-Mills coupling. Figure shows the 
overall exclusion contours for the three vector couplings. 



X. CONCLUSIONS 



We have presented 95% C.L. upper limits on the pair production of leptoquarks that 
decay to the eejj, evjj and vvjj final states. For scalar leptoquarks, the limits on the 
cross section provide lower limits on the scalar leptoquark mass of 225 GeV/c^ for (3 = 1, 
204 GeV/c^ for /3 = i, and 79 GeV/c^ for /? = 0. We have also set mass limits for vector 
leptoquarks for different couplings and have presented exclusion contours on (3 and Mlq. At 
the 95% C.L., our results exclude an interpretation of the HERA high-Q^ excess as s-channel 
scalar leptoquark production for Mlq < 200 GeV/c^ and (3 > 0.4. These results can be also 
used to set limits on the pair production of any heavy scalar particle that decays into a 
lepton and a quark as expected in a variety of models and to restrict any new leptoquark 



models containing additional fermions 
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FIG. 21: The 95% C.L. upper limits on cross sections for vector leptoquark pair production from 
all three channels combined for /3 = ^, and the LO predictions, as a function of leptoquark mass. 




FIG. 22: The 95% C.L. lower limit on the mass of first-generation scalar leptoquarks as a function 
of f3 for the individual eejj, eujj, and I'lyjj channels, and for the combined analysis. 



TABLE X: Limits on the masses of first-generation leptoquarks. 



(3 


Scalar 
(GeV/c2) 


Minimum Cross Section 
(GeV/c2) 


Minimal Vector 
(GeV/c2) 


Yang-Mills 
(GeV/c2) 


1 


225 


246 


292 


345 
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204 


233 


282 


337 





79 


144 


159 


206 



37 



1 




400 



FIG. 23: The 95% C.L. lower limit, as a function of (5, on the mass of first-generation vector 
leptoquarks with Yang- Mills couplings from the individual eejj, evjj, and vvjj channels and for 
the combined analysis. 
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